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Abstract 



The Gamow vector description of resonances is compared with the S-matrix 
and the Green function descriptions using the example of the square barrier 
potential. By imposing different boundary conditions on the time indepen- 
dent Schrodinger equation, we obtain either eigenvectors corresponding to real 
eigenvalues and the physical spectrum or eigenvectors corresponding to com- 

■ plex eigenvalues (Gamow vectors) and the resonance spectrum. We show that 

■ the poles of the S matrix are the same as the poles of the Green function and 
are the complex eigenvalues of the Schrodinger equation subject to a purely 

Oh! outgoing boundary condition. The intrinsic time asymmetry of the purely 

^ \ outgoing boundary condition is discussed. Finally, we show that the proba- 

cy ' bility of detecting the decay within a shell around the origin of the decaying 

^ • state follows an exponential law if the Gamow vector (resonance) contribution 

to this probability is the only contribution that is taken into account. 

■ 

•i-H - 

X 

H ■ I. INTRODUCTION 

Most elementary particles are only quasistable states decaying through various inter- 
actions and thus have finite lifetimes of various orders of magnitude.0 Several theoretical 
descriptions have been proposed for quasistable states. Three of the most widely used de- 
scriptions are the S matrix, the Gamow vector, and the Green function. We shall use the 
example of nonrelativistic potential scattering, in particular the square barrier potential, to 
demonstrate the connection between these three descriptions. 

Experimentally, resonances often appear as peaks in the cross section that resemble the 
well-known Breit-Wigner distribution. The Breit-Wigner distribution has two characteristic 
parameters: the energy Er at which the peak reaches its maximum, and its width Tr at 
half-maximum. The inverse of Vr is the lifetime ot the decaying =tateS The peak ot the 
cross section with Breit-Wigner shape is related to a first-order pole of the S matrix in 
the energy representation S{E) at the complex number zr = — iV {i/2. The theoretical 
expression of the cross section in terms of S{E) fits the shape of the experimental cross 
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section in the neighborhood of En (see, for example, Ref. This is why the first-order 
pole of the S matrix is often taken as the theoretical definition of a resonance. The Green 
function description treats a resonance as a pole of the Green function when analytically 
continued to the whole complex plane. 

Although a resonance has a finite lifetime, it is otherwise assigned all the properties that 
are also attributed to stable particles, such as angular momentum and charge. For example, 
a radioactive nucleus has a finite lifetime, but otherwise it possesses all the properties of 
stable nuclei; in fact, it is included in the periodic table of the elements along with the stable 
nuclei. Therefore, it seems natural to seek a theoretical description that provides "particle 
status" to quasistable states. The description of a resonance by Gamow vectors allows us 
to interpret resonances as autonomous exponentially decaying physical systems. 

The energy eigenfunction with complex eigenvalue was originally introduced by Gamow 
in his paper on a decay of atomic nuclei,@ and used by a number of authors (see, for example, 
Refs. |5|-p!T| and references therein). The real part of the complex eigenvalue is associated 
with the energy of the resonance, and the imaginary part is associated with the inverse 
of the lifetime. Gamow eigenfunctions have an exponentially decaying time evolution, in 
accordance with the exponential law observed in a decay of radioactive nuclei.&^l 

Gamow's treatment was heuristic though, and could not be made mathematically rigor- 
ous within the Hilbert space theory, because self-adjoint operators on a Hilbert space can 
only have real eigenvalues. A rigorous mathematical treatment of Gamow vectors needs 
an extension of Hilbert space to the Rigged Hilbert Space (RHS). The RHS was first 
introduced in physics in order to justify Dirac's bra-ket formalism. In RHS language, 
Gamow vectors are eigenkets of a (dual) extension of the self-adjoint Hamiltonian. This 
extension can surely have complex eigenvalues. Using the RHS formalism, one can prove 
that the time evolution of the Gamow vectors is governed by a semigroup, expressing time 
asymmetry on the microscopic level. One can also obtain an exact golden ru iJ for the decay 
of the state described by the Gamow vectors. This golden rule reduces to the Fermi-Dirac 
golden rule in the Born approximation. 

There are several pedagogicalpapers on Gamow vectorSii and related topics such as a 
decay,i'@ barrier penetration,!!'!! and exponential decay!!^ (this list is not exhaustive). 
Some of these papers consider alternatives to the Gamow, S matrix, and Green function 
descriptions. Especially suggestive is the approach of Refs. ^ and ^ where among other 
methods Feynman's path integral is used. In this paper, we discuss at the graduate student 
level some of the issues not covered by Refs. 



In Sec. [II A| , we calculate the Dirac kets for the square barrier potential. These Dirac 
kets are monoenergetic eigensolutions of the time independent Schrodinger equation. They 
are not square normalizable, and therefore they cannot represent a wave packet. Actually, 
they are members of a continuous basis that expands the wave functions in the Dirac basis 
vector expansion. 

The well-known expression for the S matrix in the energy representation is provided in 
Sec. 11 Ij . The resonances will be defined as the poles of the S matrix. In Sec. [11 (J| , we 



calculate the Gamow vectors as the solutions of the time independent Schrodinger equation 
with complex eigenvalues subject to purely outgoing boundary conditions. 

The Green function and its poles are calculated in Sec. [II D[ It will be apparent that the 
poles of the Green function are the same as the poles of the S matrix and are the complex 
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eigenvalues obtained from the purely outgoing boundary condition. The residue of the Green 
function at the resonance energy is expressed in terms of the Gamow eigenfunction. 

In Sec. |11E| , the basis vector expansion generated by the Gamow vectors is provided. 
Sections [III A| and PIB| discuss the time asymmetry built into the purely outgoing boundary 
condition. Section [111 C| treats the exponential decay law of the Gamow vectors. 

In the paper, there will be some mathematical interludes enclosed by brackets. These in- 
terludes are not essential to understand the paper (although they are essential to prove other 
fundamental resultsSlli) and may be skipped in the first reading. They are recommended 
to readers who are seriously interested in the mathematical framework that supports the 
Gamow vector approach to resonance scattering. 



II. RESONANCES FOR A SQUARE BARRIER POTENTIAL 

A. Time Independent Schrodinger Equation 

We consider a three-dimensional square barrier potential of height Vq > and calculate 
the energy kets using the time independent Schrodinger equation 

H\E)=E\E). (1) 

Equation (0) is an eigenvalue equation of the Hamiltonian H. When the eigenvalues E 
belong to the continuous spectrum of H, then the solutions \E), also called Dirac kets, are 
given by vectors that lie outside the Hilbert space, that is, they are not square integrable. 
Following von Neumann, these kets are sometimes interpreted as states that are "very near 
a proper state. "il Following Dirac,0 we interpret these kets as members of a continuous 
basis system that spans the space of wave functions. Therefore, every wave function if can 
be expanded in terms of the Dirac kets \E, I, m) = \E) as 

^ ^ POO 

^ = J2J2 dE\E,l,m){EJ,m\ip). (2) 
1=0 m=-r^ 

[Although Dirac introduced Eq. (0) on heuristic grounds, its mathematical rigor was later 
established by the Gelfand-Maurin theoremllll (also called the nuclear spectral theorem) 
within the Rigged Hilbert Space (see also Refs. p!6|JT8| - pT| , p2| , and | 



To calculate the possible set of (real) eigenvalues and their corresponding eigenvectors 
in our example, we solve Eq. (|T|) in the position representation, 

{x\H\E) = V + Vix)^ {x\E) = E{x\E) , (3) 

where is the three-dimensional Laplacian and 

<r < a 

V{x) = V{r) ={Voa<r<b (4) 
6 < r < oo . 
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Because V{x) is spherically symmetric, we can use spherical coordinates x = (r, 6, 0) to 
solve Eq. (|^), which in spherical coordinates reads 

(r, /, m) = [^-g^r + + V^W) (r, ^, /, m) = E{r, 6, /, m) . 

(5) 

By separating the radial and angular dependences, 

(r, 9, ^\E, I, m) = {r\E)i {9, </)|/, m) = -xi{r- E)Yi,^{9, </>), (6) 

r 

where Yi^^ni^^, 4') are the spherical harmonics, we obtain for the radial part 

(^^, + ^^f^ + Vir))xdr; E) = Ex^r; E) . (7) 

In this section, we shall restrict ourselves to the case of zero orbital angular momentum (the 
higher-order case is treated in the appendix). We write X/=o('"; E) = x(r; E) and obtain 

- ^-^Xir; E) + Vir)xir; E) = Ex{r; E) . (8) 

The solutions of Eq. for each value of E are the kets (more precisely, the radial part of 
the energy kets in the position representation {x\E) = (r, 9, (p\E, 0, 0) = ^^^^^^ Yq^o(9, 0) ). In 
our case, 

( aie'^'' + /?ie-*'' < r < a 
X(r; E) = I a2e'^' + p2e-'^' a<r <h (9) 
\ J^ie'^"- + J^2e-'^'' b<r <oo, 

where 



k=\l'-^E (10) 



is the wave number of the particle, and 



Q=^k^--'i^V„=f-^(E-V„). (11) 

The coefficients a, f3, and are functions of k (and therefore of the energy E). Their 
possible values are restricted by the boundary conditions that we need to impose on the 
solutions (H) of the radial Schrodinger equation. Some of these boundary conditions are due 
to the fact that at least the second derivative of x must be well defined, which is why x 
must be a continuous function of r with a continuous derivative. Another boundary condition 
makes the solution vanish at the origin (see Eq. ( |12a|) ). This condition can be understood 
by viewing the analogous one-dimensional problem of Eq. (^. In the one-dimensional case 
the potential is infinite for r < 0. The function T2e~^^^ represents an incoming probability 
wave of amplitude and T\t^^'^ represents an outgoing wave of amplitude T\. Because 
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there cannot be any transmission into the region r < 0, x must vanish at the origin.a 
Mathematically, this condition is related to the self-adjointness of the Hamiltonian.iilli 
Finally, the eigenfunction x{^', E) is assumed to be bounded (see Eq. (|121|)). We impose this 
boundedness condition because we want the set of eigenvalues to coincide with the Hilbert 
space (that is, the physical) spectrum of the Hamiltonian. In our case, the Hilbert space 
spectrum is [0, cxd) (see Refs. ^ and |16D. Then, the boundary conditions that we impose 
upon Eq. (13) read 



X(0;E) = 
x{a-0-E) = x{a + 0;E) 
x'ia-0;E) = x'ia + 0;E) 
x{b-0;E) = xib + 0;E) 
x'{b-0-E) = x'{b + 0-E) 
\x{r;E)\ < oo. 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 
(12f) 



If E' is a negative real number (or complex), then there are no coefficients of x('"; E) in (^ 
for which the boundary conditions ([T^ can be satisfied unless they are trivially zero. To be 
more precise, if E' is a complex or a negative real number, the corresponding eigensolution 
x{r', E) of Eq. (^ does not satisfy the boundary condition (|12fl) , even though it satisfies the 
other boundary conditions ( |12a| )-( [T^ ). If is a real number in [0,cxd), the corresponding 
eigenfunction satisfies all the boundary conditions in (|12[). Therefore, the boundedness of 
the eigenkets is what forbids the negative (and the complex) energies and singles out the 
physical (Hilbert space) spectrum. 

If we use the notation a = 2iai, we rewrite the boundary conditions ([l^) as the following 
equations for the coefficients: 



aae^^'^ + Z^se-'^'^ = asin(A;a) 
iQ{a2e'^'' - p2e~'^'') = akcos{ka) 

ik^Tie''''' - J^ae"''^'') = iQ{a2e''^'' - (326-'^') 
After straightforward but tedious calculations, we find that 

k 1 
sin(A;a) H — — ■ cos(A;a) a(k) 
iQ J 

k 1 

sm(ka) cos(A;a) a{k) 

iQ J 

k 



a2{k) = 


-'- ~iQa 

2 


m) = 


2 




^—ikb - 


^l(^) = 


4 [ 



sin(/ca) + — cos(A;a)) 
iQ 



+ (1 - ^)e~*«('-'^)(sin(A;a) - -^cos(A;a)) aik) 
k iQ 



J'2{k) 



k' 



sin(A;a) + — cos(A;a)) 



tQ 

(1 + ^)e-'^'^''~^\sm(ka) - cos( ka))] aik) . 
k iQ J 



(13a) 
(13b) 
(13c) 
(13d) 



(14a) 
(14b) 

(14c) 

(14d) 
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Thus, for each E in [0, oo), there exists a solution of the eigenvalue equation in the position 
representation 



(r, 6, m = Fo,o(^, 0) = a/^ , < E < oo (15) 

r r V 47r 

with x(r; E) given by 

{a(fc) sin(/cr) < r < a 

a2{k)e'^'' + /?2(A;)e-^^^ a < r < 6 (16) 

J'i(A;)e^'=^ + J'2(A;)e-^'=^ 6 < r < oo . 

The coefficients a2(fc), P2{k), and J-'i^2{k) in Eq. (|16D are given by Eq. (p^) in terms of a{k). 
Usually, a{k) is chosen such that the kets in Eq. ( [16|) are 5-normalized, although in the 
following sections we shall assume that a{k) = 1. 

Equation (^) means that for each energy eigenvalue E in the spectrum [0, oo) of the 
Hamiltonian, there exists an eigensolution of the Hamiltonian satisfying the boundary con- 
ditions ([T2|) . We can expand any wave function (p in terms of these eigenf unctions (for / = 0) 
as 

^(r, e, 0) = rdE^^^^Yo,oie, 0) ifiE) , (17) 
Jo r 

or in bra-ket notation 

/•oo 

(r,^,0|y.) = / dE{r,9,<P\E){E\v)- (18) 
Jo 

We can interpret (|l^) by saying that any wave function (y9 is a continuous linear superposition 
of the eigenfunctions (r, 6,(j)\E) in Eq. (^). Each eigenfunction is weighted by {E\(f), which 
represents the wave function in the energy representation. 

[Mathematically, the eigenkets {r,6,(f)\E) allow us to go from the energy representation 
(p{E) = {E\ip) to the position representation ip{x) = {r,9,(j)\ip) and vice versa, that is, 
they are continuous transition matrix elements. Because the monoenergetic {r,6,(f)\E) are 
not square integrable, they are not in the Hilbert space. Thus Hilbert space methods are 
not sufficient to handle them, and an extension of those methods is needed. As shown in 



Refs. 16,n8lp2l and 33, the extension that seems to be the most convenient is the Rigged 



Hilbert Space. 



B. S-Matrix Approach 

We consider now the scattering process of a particle beam off the square barrier poten- 
tial. The S matrix for this process relates the incoming (or prepared) and the outgoing (or 
detected) beams. We write the S matrix in the energy and angular momentum representa- 
tion, 

{-E, /, m\E\ /', m'+) = SiiE) 5{E - E') 5i,, 6m,m' , (19) 
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where \E,l,m^) are the kets that solve the Lippmann-Schwinger equation, and {^E,l,m\ 
are their corresponding bras. In the position representation, these Lippmann-Schwinger 
eigenfunctions are given by (see Refs. 



{r\E^) = Zi%^ (20a) 
{r\E-) = 1^^. (20b) 
The probabihty amphtude of detecting an out-state ip~ in an in-state is (see Refs. ^ and 

m 

{^-,^+) = J2J2 dE{^p-\E,l,m-)SiiE){+E,l,m\^+). (21) 

/=0 m=-l 

In this section, we restrict ourselves to the case of zero angular momentum, 1 = 0, that is, 
to the first term of Eq. (pT]). The partial S matrix for the case of zero angular momentum 
will be denoted by S{E) = Si=o{E). 

Because the S matrix relates incoming and outgoing waves far outside the interaction 
region, it suffices to focus on the region r > b. In this region, we have an incoming spherical 
wave e~'^''/r with amplitude J-'2{k) and an outgoing spherical wave e*'^'"/r with amplitude 
J-'i{k). Therefore, the S matrix in the energy representationSH is given by the ratio 

where k is given by Eq. ([l^). It is easy to verify that the S matrix is unity when there is no 
potential, 

lim S(k) = 1 , (23) 

and that S*{k)S{k) = 1, that is, |5'(A;)| = 1 for every k. 

The S matrix S{k) in Eq. ( [2^ ) is well defined for each real value of k, because its 
denominator is never zero when k is real. The / = resonances are associated with the poles 
of the analytic continuation of S{k) into the entire complex plane. The relation k'^ = 2mE/h'^ 
provides a Riemann surface for this extension in a natural way. The analytic continuation of 
the numerator and the denominator of S{k) yield two analytic functions J-'i^2{k)- Therefore, 
the continuation of S{k) is analytic except at its poles. These are precisely the zeros of the 
denominator of S{k) (see Ref. 



^2(fc) = 0, (24) 

which leads to 

(1 _ ^)e^Q(fe--) [sin(A;a) + 4 cos{ka)] + (1 + 9.)e-'<^ib-) [sin{ka) - A cos(A;a)] = . (25) 
k iQ k iQ 

The solutions of Eq. (|25|) are the (S'-matrix) resonances of the square barrier potential. 
Equation (|25|) has a denumerable infinite number of complex resonance energy solutions. 
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These solutions come in pairs ER±iTR/2 (see Fig. |l]). The pole ER—iVR/2 is associated with 
the decaying part of the resonance, and is located on the lower half-plane of the second sheet 
of the two-sheeted Riemann surface corresponding to the square root mapping (see Fig. |I]a). 
The pole Eji + iTji/2 is associated with the growing part of the resonance, and is located on 
the upper half-plane of the second sheet of the Riemann surface (see Fig. ^). In the wave 
number plane, this pair of energy poles corresponds to a pair of poles ±Re(/c) — zlm(fc) in 
the lower half-plane that are mirror images of one another with respect to the imaginary 
axis (see Fig. ||). 

The width of the resonance increases as the energy increases, and therefore the lifetime 
r = l/F/j decreases. The resonances whose energies are below the threshold E = Vq are 
close to the real axis. As E increases, the resonances move away from the real axis toward 
infinity. 



C. Gamow Vector Approach 

In this section, we determine the state vectors of the resonances of the square barrier 
potential (|) as solutions of the time independent Schrodinger equation with complex eigen- 
values. The state vectors will be eigenvectors as in Eq. (|I]). Because in Sec. P A| we found 
all the eigenvectors of H with the boundary conditions ([l2|) and determined that they all 
have positive real eigenvalues, we will have to change these boundary conditions to purely 
outgoing boundary conditions.!! We will see that there is a one-to-one correspondence be- 
tween the complex poles of the analytically continued S matrix and the complex eigenvalues 
obtained from purely outgoing boundary conditions. These eigenvectors of the Hamiltonian 
with complex eigenvalues, which are associated with the poles of the S matrix, are called 
Gamow vectors or Gamow kets. 

The time independent Schrodinger equation for the Gamow vectors is given by 

H\zi,) = zr\zr) , (26) 

where zr and \zr) represent the complex resonance eigenvalue and the Gamow vector, 
respectively. In the position representation, Eq. (HR) becomes 



or 



{x\H\zr) = zr{x\zr), (27) 



+ Vix)^ {x\zr) = zr{x\zr) . (28) 



momentum, we obtain 



After writing Eq. (|28|) in spherical coordinates and considering only the case of zero angular 

2m rfr2-^*^^' ^^"^ ^ ^(^)x('"; zr) = ZRxir; zr) . (29) 

The solution of Eq. (p9D for each value of zr is the s-Gamow vector (more precisely, the 
radial part of the s-Gamow vector in the position representation {x\zr) = (r, 6, (f)\zR, 0, 0) = 
1/r x(r; ZR)Yofi{9, 0)). The solution of Eq. (^) has the same form as Eq. (H) 
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aie*'"' + Pie''-'''' < r < a 
X{r; Zr) = { 026'^'' + p2e-'^'' a<r <b (30) 



but now the wave number 



of the resonance is complex, and 



/2m 

k=\I^ZR (31) 



Q = \jk--^-^yo = f-^{zR-Vo). (32) 

Because the second derivative of xij] zr) must be well defined, the solution ( |30D must be 
continuous with continuous derivatives as in (|12bD -( [T2e|) . For the Gamow vector we shall 
also choose the condition that the radial part xl?"! zr) vanishes at the origin as in Eq. (|12aD . 



However, at infinity we choose the purely outgoing boundary condition, which means that 
far from the potential region, the solution reduces to an exponential of the type e*'^^, but 
not of the type e~**"'. If k is complex, this purely outgoing boundary condition does not 
mean that there are only outgoing waves. In fact, we have outgoing waves only when Re(/c) 
is positive, and incoming waves when Re(/c) is negative. Thus the boundary conditions that 
we impose on the Gamow vectors are 

(33a) 
(33b) 
(33c) 
(33d) 
(33e) 

^ 00 . (33f ) 

The purely outgoing boundary condition (|331|) is often written as 





X(0; Zr) 


= 




-0;zr) 


= xia + 0; Zr) 


X{a 


- 0; zr) 


= x'ia + 0; Zr) 


Xib 


- 0; Zr) 


= x{b + 0; Zr) 


x'{b 


- 0; zr) 


= x'ib + 0;zR) 




Xir; Zr) 


Akr 

~ e , 



lim MjiM _ ^kx{r■ zr) = . (34) 



+00 



One can easily check that Eq. (|3^ is equivalent to Eq. ( p3f|) . 

The boundary condition ( |331]) leads by Eq. (|30|) to T2 = 0. Because this condition is the 
same as the condition (^) for the complex poles of the S matrix (122|) , the set of complex 
eigenvalues zr must include the set of complex solutions of Eq. ( pSj) (which are the S'-matrix 
resonance poles). We shall show that these two sets of solutions are the same. 

If we define a = 2iai, the boundary conditions (^) can be written in terms of the 
coefficients as 

ase^^" + ^se"^'^" = a sin(fca) (35a) 
iQ{a2e'^'' - /^se"^^") = ak cos{ka) (35b) 
J^,e'^^ = a2e'^^ + /^se"'^^ (35c) 
ikj^^e'''^ = iQ{a2e'^^ - p2e-'^^) ■ (35d) 
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If we write this set of linear equations as a matrix equation, we obtain 



/ sin(A;a) 
k cos{ka) 


\ ike'''^ 








-iQe'^" iQe 



oiQb 



-iQb 
iQb 



/0\ 






(36) 



This is a homogeneous system of four equations with four unknowns. The system has a 
non-trivial solution if, and only if, the determinant of the coefficients is equal to zero. 



sin(/ca) 





gjQa 


— 


-iQa 


k cos{ka) 







iQe~ 


-iQa 





^ikb 


_^iQb 


—e" 


-iQb 





ike^^^ 




iQe 


'iQb 



0. 



Straightforward calculations lead to 
(1 - ^)e'«('-'^)[sin(fca) + ^ cos(A;a)] + (1 + ^)e-'^('-'^) [sin(A;a) - A cos(A;a)] 



(37) 



0. (38) 



Thus we obtain exactly the same condition as the condition ( P5|) for the poles of the 5* 
matrix. If we compare the boundary conditions (|T2p imposed upon the Dirac kets with the 
boundary conditions (|33D imposed upon the Gamow vectors, we can see that the boundary 
condition that singles out the resonance spectrum is the purely outgoing boundary condition 

As we mentioned earlier, the solutions of Eq. (|38|) (or Eq. (^)) always come in pairs. 
The eigenvector corresponding to the complex eigenvalue E^i — iT ji/2 is the decaying Gamow 
vector in the position representation, whose radial part, up to a normalization factor, is 



sin(/crfr) < r < a 

a2{kd)e'^'^'' + p2{kd)e-'^'"' a<r <b 
J^i{kd)e'^'^'' h <r < oo 



(39) 



where kd = ^J2mir9 [ER-iYnjl) and Q\ = k\- 2m/h^ Vq. 



For Eji + iTji/2, we obtain the growing Gamow vector in the position representation. Its 
radial part, up to a normalization factor, is 



X^™s(r;4) 



sm{kgr) 
a2{kg 



JQnr 



Ti{kg)e'^^'- 



+ P2ikg 



< r < a 
a < r < b 
b < r < oo 



(40) 



where kg = ^^mlfi^ {Er + iTR/2) and Q^ = - 2m/ Vq. 

[Gamow eigenfunctions are not square integrable, that is, they do not belong to the 
Hilbert space. Thus, as Dirac kets, they must be handled by methods that extend those 
available in the Hilbert space framework. The Rigged Hilbert Space provides those methods 
(see for example, Refs. and 18).] 
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D. Green Function Approach 



In this section, we perform the calculations for the resonance energies using the Green 
function method. The radial Green function satisfies (see for instance, Refs. ^ and ^) 

+ - G(r, r'; E) = -6{r - r% (41) 

subject to various boundary conditions. Equation (|4l| ) says that for r 7^ r', G{r, r'; E) obeys 
the radial Schrodinger equation (^). At r = r' it is continuous, but its derivative has a 
discontinuity due to the delta function, 

|-G(r' + 0, r'- E) - |-G(r' - 0, r'; = . (42) 

We want G to satisfy the same kind of boundary conditions at zero and at infinity that we 
imposed on the other methods. We shall also require G to be regular at zero and a Durely 
outgoing boundary condition at infinity. The expression for this Green function isii'0 

Gir r'- E) = ^ ^(loWlr^ 

where r< refers to the minimum and r> to the maximum of r and r'. The function x is 
the solution of Eq. (§) that vanishes at the origin, ip is the solution of Eq. (|^) that satisfies 
a purely outgoing boundary condition at infinity, and W{x,ip) is their Wronskian. In this 
case X is just the solution in Eq. (p^Bf) . The function ip satisfies Eq. (|]) and the boundary 
conditions 

^{a - 0; E) = ipia + 0; E) (44a) 
^'{a -0]E) = tp'{a + 0; E) (44b) 
^{b - 0; E) = tPib + 0; E) (44c) 
ij'ib -0;E)= ^p'ib + 0; E) (44d) 
ij{r;E) e'^\ r^oo, (44e) 

The solution ip is then given by 

{ai(/c)e*'"' + 6i(A;)e"*'"' 0<r<a 

a2(fc)e*'3^ + 62(A;)e-^«^ a<r <b (45) 

e*'^'' 6 < r < 00 , 

where the coefficients a{k) and b{k) are functions of k that make ip satisfy the boundary 
conditions of Eq. (^^ . The Wronskian of x S'Ud ip is 



W{x.^){k) = x{rW{r) - x'(r)^(r) = 2ikJ'2{k) . (46) 
From Eqs. ( P5| ) and (^), we obtain the Green function, 

Gir r'- k) = ^ Xir<;k)^piry,k) 
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The poles of this Green function are simply the zeros of its denominator.E§ Those poles are 
the same as the S*- matrix poles (^51) and as the complex eigenvalues in Eq. (^). 

It is worth noting that the Green function (^) is an outgoing Green function only when 
Re(/c) is positive. If Re(/c) is negative, the Green function (^7]) has a purely incoming 
character (cf. Sec. |T|). Also, as k approaches the real positive k axis or as E approaches the 
right-hand cut from above, the Green function G becomes the outgoing Green function G"*"; 
as k approaches the negative real axis (from above) or E the right-hand cut from below, G 
becomes the incoming Green function G".!! 

The residue of the Green function (^7p at a decaying resonance wave number is given by 

2m 1 

|G(. r'; n].... = ___^_^,-.».(.,; . (48) 

A similar relation between the residue of G(r, r'; k) at the growing wave number pole kg and 
the growing Gamow eigenfunction (P0[) holds. 



E. Complex Basis Vector Expansion 

The Dirac kets are basis vectors that were used to expand a wave function Lp as in 
Eq. (|l^). This expansion is known as the Dirac basis vector expansion.^ The Gamow 
vectors are also basis vectors. The expansion generated by the Gamow vectors is called the 
complex basis vector expansion. However, the Gamow vectors do not form a complete basis. 
The complex basis vector expansion needs an additional set of Dirac kets corresponding to 
the energies that lie in the negative real axis of the second sheet of the Riemann surface.^" 
This has been realized by other authorJ30'0 who used the Green function approach. In 
this section, we shall expand a wave function in terms of the Gamow vectors (which contain 
the resonance contribution) and a continuous set of Dirac kets (which is interpreted as a 
background contribution). 

[Rigorously speaking, the complex basis vector expansion is not valid for every normaliz- 
able wave function, that is, for every element of the Hilbert space, but only for those square 
integrable functions that are also Hardy functions. 

Experimentally, we can only measure the probability for an event to take place. In a 
scattering experiment this is the transition probability from an in-state ip'^ to an out-state 
ip~ . The amplitude of this probability is calculated from the scalar product between 
and ^Z'". In terms of the quantities that we already have calculated, this amplitude is given 
byiS 

/■oo 

(7/.-,^+) = / {r\E~)S{E){E+\^+)dE. (49) 
Jo 

We now extract the resonance contribution from Eq. (|i9|). This resonance contribution is 
carried by the Gamow vectors. In order to do so, we deform the contour of integration into 
the lower half-plane of the second sheet of the Riemann surface for the S matrix, where the 
resonances poles are located (see Fig. ^). If we use the results in Ref. we can write 
Eq. (IID as 



12 



•^0 n=0 

where z^^n = En — iTn/'2. denotes the nth decaying pole and r„ denotes the residue of S{E) 
at Zd^n- In Eq. (0), the integration is done on the negative real semiaxis of the second sheet 
of the Riemann surface. The series in Eq. (|50D can be shown to be convergent.@ Omitting 
in (|50|), we obtain the complex basis vector expansion for the in-states, 

/ \E+){+E\^+)dE - 27rz^r„|^^-„)(+^d,„|(^+) . (51) 

Jo n=0 



In Eq. (plD, the infinite sum contains the resonances contribution, while the integral is 
interpreted as the background contribution. 

Similarly, we obtain the complex basis vector expansion for the out-state but now 

we deform the contour of integration into the upper half-plane of the second sheet of the 
Riemann surface, where the growing resonance poles are located (see Fig. ^d). The result is 



/• — CO '■^ 

/ \E-){E- \r)dE + 27:1 J2 {-^;,n\r) , (52) 

■^0 „=o 



where z*^^ = En + «r„/2 is the nth growing pole, and r* is the residue of S{E) at 2* „. The 
integration in Eq. (|5^) is performed on the negative real semiaxis of the second sheet of the 
Riemann surface. The series in Eq. (B^) has been shown to be convergent .0 



III. TIME ASYMMETRY OF THE PURELY OUTGOING BOUNDARY 

CONDITION 

The semigroup time evolution of the Gamow vectors expresses the time asymmetry built 
into them.iHil0 We will show here that the purely outgoing boundary condition that singles 
out the resonance energies also has an intrinsic time asymmetry. To be more precise, we 
will show that the purely outgoing boundary condition should be read as purely outgoing 
only for the decaying part of a resonance and as purely incoming for the growing part of 
the resonance. Obviously, the purely incoming boundary condition is the time reversed of 
the purely outgoing one. Therefore the growing Gamow vector can be viewed as the time 
reversed of the decaying Gamow vector. 



A. Outgoing Boundary Condition in Phase 

First, we study the meaning of the purely outgoing boundary condition when it is imposed 
on the decaying part of the resonance. The complex energy associated with the decaying 
part of a resonance is Zd = Eji — iTfi/2 {Eji, F/j > 0), which lies in the lower half-plane of 
the second sheet of the Riemann surface (see Fig. |I]a). Its wave number kd = Re(fc) — zlm(/c) 
(Re(fc), Im(fc) > 0) lies in the fourth quadrant of the wave number plane (see Fig. 
The decaying Gamow vector T^^ecaymg ( |gg[ ) was obtained after imposing the purely 
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outgoing boundary condition ( P31]) on Eq. (|30D. If we had not imposed this condition, we 
would had obtained a solution of the form (|30|) , and every complex number would had been 
an eigenvalue of the Hamiltonian. In the region r > b, this solution would had been the sum 
of two linearly independent solutions: 

decaying f\ = JT ^-ik^r -iz^t/h 
Aincoming V 5 ''7 2<- 

= (jFse"^"^''^''"'^^*/^^'^^) e-*^'^^'')^-*-^^*/^ , r > b, (53) 
which we call the incoming decaying Gamow vector, and 

decaying/ ,n _ r ik^r -iz^t/h 
Aoutgoing V ) 1^ ^ 

= (^JZ^(^Hk)r~VRt/{2h)^^ ^V.e{k)r-iERt/h ^ r>6, (54) 

which we call the outgoing decaying Gamow vector. These names come from the standard 
interpretation of plane waves with a complex exponent (see for instance Ref. ^): the expo- 
nential with a purely imaginary exponent (the term that carries the phase) is interpreted as 
the term that governs the propagation of the wave, and the exponential with the real expo- 
nent is interpreted as the term that just changes the amplitude of the wave on the surfaces 
of equal phaseB We are going to interpret Eqs. ( ^3|) and (54) in the same fashion. The 



terms between parentheses in Eqs. (^3|) and (|5^ ) determine the amplitude of the waves. The 
propagation of Xoutgoing governed by ^'^'^'^(ky^-tEiit/h^ therefore Xoutgolng an outgoing 
wave (in phase). Analogously, the propagation of xfncoming governed by e-«Re(A:)rg-j_BBt/ft^ 
and thus xfncoming is incoming wave (in phase). Imposing the purely outgoing boundary 
condition JF2 = is tantamount to forbidding xfncomijig- Thus for the decaying part of the 
resonance, the purely outgoing boundary condition allows only purely outgoing waves. 

The meaning of the purely outgoing boundary condition applied to the growing part of 
the resonance is the opposite. The growing energy eigenvalue Zg = E^ + iV r/2 [Er, > 0) 
lies in the upper half-plane of the second sheet of the Riemann surface (see Fig. ^), and its 
wave number kg = — Re(A;) — zlm(/c) (Re(/c), Im(fc) > 0) lies in the third quadrant of the wave 
number plane (see Fig. Q). The growing Gamow vector in Eq. (^Dj) was obtained 

after imposing the condition ^2 = on Eq. (|30D. If we had not imposed this condition, in 
the region r > b, the solution would had been the sum of two linearly independent solutions: 

growing / ,^ ^ -p^ikg'r -izgt/h 
AincomingV ) 1^ ^ 

= (jrie^™(''^'"+'^^*/^^^)) e-*^«=(^)''-*'^«*/^ , r > b, (55) 
which we call the incoming growing Gamow vector, and 

growing / ,^ _ -ik^r -izgt/h 

AoutgoingV ' ~ 2C c 

= (jr2e"^™('')''+'^«*/(2'^)) (^'^<k)r-iERt/h ^ r > b, (56) 

which we call the outgoing growing Gamow vector. The names also come from the standard 
interpretationl^i of plane waves with a complex exponent. Therefore, the purely outgoing 
boundary condition JF2 = 0, when applied to the growing part of a resonance, bans Xoutgoing 
and allows only purely incoming waves. 
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B. Outgoing Boundary Condition in Probability Density 



In Sec. [IIl7\| , we showed how the time asymmetry built into the purely outgoing boundary 
condition affected the phase of the Gamow vectors. In this section, we show the same time 
asymmetry but now consider the probability density of the Gamow vectors. 

For the decaying part of the resonance, the probability densities (before imposing the 
purely outgoing boundary condition) are obtained by taking the square of the absolute value 
of Eq. 

decaying/ ,n ^ i decaying / ,n|2 ^ i -77 1 2 -2Im(fe)r-rflt/^ 
rincomingV 5 I Aincoming V ' I K2| ^ 

= |JSre-r^/^(*+^/'') , r>6, (57) 

which we call the incoming decaying probability density. Similarly, the square of the absolute 
value of Eq. (|5^) yields 



decaying/ ,n ^ 1 decaying/ f\\2 ^ \ -p \2 2lxii{k)r ~V / h 
routgoingV ' I A outgoing V 5 I Kl| ^ 



|^^|2g-r«A(t-r/.)^ r>6, (58) 



which we call the outgoing decaying probability density {v = r^/(2Mm(fc))). By imposing 
the purely outgoing boundary condition J^2 = 0, we allow only Eq. (|5^) and forbid Eq. {^7\j, 
which we interpret by saying that we have a purely outgoing probability density condition 
for the decaying part of the resonance. 

For the growing part of the resonance, the probability densities (before imposing JF2 = 0) 
are the square of the absolute value of Eq. (^) , 

growing / ,^ ^ i growing / ,^,2 ^ \p \2 21m{k)r+rRt/h 
rincomingV ' ''Z I AincomingV ' I I'' 1| '-' 

= |J-i|V«/^(*+'-/'^), r>b, (59) 
which we call the incoming growing probability density. Similarly, from Eq. (0), we have 

growing / ,\ _ I growing / ,\|2 _ 1^7 \2 -2lm{k)r+rRt/h 
^outgoing V ' I AoutgoingV ' I r 2| c 

= |J^2| , r>b, (60) 

which we call the outgoing growing probability density. For this growing part, the condition 
JF2 = only allows waves with purely incoming probability densities. 

In short, the purely outgoing boundary condition ( ^3iD must be read as purely outgoing 
(in phase or in probability density) only for the decaying part of the resonance and as purely 
incoming (in phase or in probability density) for the growing part of the resonance. 



C. Exponential Decay Law of the Gamow Vectors 

We want to determine the probability VAroit) of detecting the decaying state within 
a shell of width Aro outside the potential region (r > b). This is the probability that 
is measured by the counting rate of a detector placed, for example, outside a radioactive 
nucleus from which an a-particle is emitted. We assume that the detector surrounds the 
nucleus completely and that is at a distance > b from the center r = 0. 
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Theoretically, the probability VAroit) to observe an in-state (p^ at time t within the 
interval Aro around the surface r = rg is given by 



Experimentally, the probability of finding the decaying state particle around ro, that is, 
the counting rate of the detector, is not defined for all times t: a resonance must be first 
prepared before the system can decay. The time at which the preparation of the resonance 
is finished and at which the decay starts can be chosen arbitrarily (we choose it to be 0). For 
example, the a particle emitted by an a-unstable nucleus travels at speed v = Tr/ {2hlm{k)) 
and reaches the point at the time t(/'o) = ro/v. For times less than t(ro), the a particle 
is not there yet, and therefore the counting rate measured by a detector placed at tq is zero 
for times t < tq/v. Whatever would have been counted by the detector before the instant 
t(ro) at To cannot be connected with the decaying state. Thus the theoretical probability 
to detect a resonance at tq should be zero for t < tq/v. This is an instance of the time 
asymmetry built into a decaying process. 

Experimentally as well, the decay of unstable systems usually follows the exponential 
law (cf. Refs. HHIl). 

[Hilbert space cannot accommodate either the time asymmetry of PAro(^)0 or the ex- 
ponential decay law.0 To account for these two features, we should use the Rigged Hilbert 
Space.iBlii In this formulation, the Gamow vectors have an asymmetric time evolution 
given by a semigroup, which accounts for the time asymmetry of a resonant process. The 
behavior of the semigroup evolution is in contrast to the time-symmetric Hilbert space time 
evolution, which is given by a group.] 

We are going to show that the exponential decay law holds if we consider only the 
resonance (Gamow vector) contribution to the probability (|6TD. In Sec. |1I E| , we used the 
Gamow vectors as basis vectors to expand the in-state in terms of the background and 
the resonance contribution (see Eq. ([5T|)). To calculate the resonance contribution to the 
probability (|6TD , we approximate by the Gamow vector by neglecting the background 
term in Eq. ([51|), 




(61) 



decaying 



(^+(r,^,0)~^^(r,^,0) 



>o,o(^,0). 



(62) 



r 



Thus the resonance contribution to the probability is 




(63) 



The time evolution of the Gamow vector &m is given by 




(64) 



and therefore 



-v' 

■iiEii-irji/2)t/h X_ 



decaying 



{r,e,<p\i,^{t))=e 



r 



(65) 
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If we substitute Eq. (^) into Eq. (|63|), we obtain 

JAro 

J Aro 

/■-ro+Aro 

= e-r«*/^|^i(A;)|2 / dre^M^K 

p2Im(fc)Aro _ 1 

_ -rRt/h\^ (i.\\2 21mik)ro l_ ^ 

2Im(fc) 

= I J-i(A;)|2Aro e-r«/'^(*-^oH , t > r^jv , (66) 
where we have used the approximation Atq small in the next to the last step. Therefore, 
VAroit) ~ |^i(A;)|2Aroe-r«/"(*-'^«/'^) , t > ro/v . (67) 



Equation ([67| ) represents the resonance contribution to the counting rate measured by a 
detector placed at tq. This resonance contribution reaches its maximum at t = tq/v and 
decreases exponentially as time goes on. Therefore, the Gamow vector (resonance) contri- 
bution to the probability VAroit) follows the exponential decay law. 



IV. CONCLUSIONS 

We have studied the different kinds of boundary conditions that we need to impose on 
the time independent Schrodinger equation in order to obtain either Dirac kets (scattering 
states) and the physical spectrum or Gamow vectors (resonance states) and the resonance 
spectrum. By imposing the boundary condition of boundedness (|121|), we obtain the Dirac 
kets and the spectrum [0,oo). And by imposing the purely outgoing boundary condition 
( P31| ) , we obtain the Gamow vectors and the resonance spectrum of Fig. |1]. This purely 
outgoing boundary condition produces the same resonance spectrum as the S matrix and 
the Green function. 

The Gamow vectors have been used as basis vectors in the complex basis vector expan- 
sions ( pTD and (|5^. However, they do not form a complete basis, and therefore a continuous 
set of Dirac kets was added to complete them. The expansions ( |51| ) and (|5^ ) extract the 
resonance contribution out of the normalized in- and out-states, respectively. 

We have uncovered the time asymmetry that arises from the purely outgoing boundary 
condition. We have seen that the purely outgoing boundary condition should be read as 
purely outgoing only for the decaying part of the resonance, and as purely incoming for the 
growing part of the resonance. 

The exponential law has been shown to hold if the background term of the complex 
basis vector expansion is neglected — only the resonance (Gamow vector) contribution to the 
probability is taken into consideration. These conclusions are not restricted to the example 
given here and are in fact applicable to a large class of potentials that includes potentials 
with compact support. 
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APPENDIX A: HIGHER ORDER ANGULAR MOMENTUM 

We extend the analysis done in Sec. ^to the case in which the angular momentum / > 0. 



1. Eigenkets 

The radial part Xi{''^i E) of (r, ^, (j)\E^ /, m) satisfies Eq. (|^. The solution of this equation 
can be written as 

(anhjikr) + I3iih~[{kr) <r < a 
anhJiQr) + ^nkiQr) a<r<b (Al) 
Tiihl{kr) + Ti2hY{kr) b < r < oc , 

where k and Q are given by Eqs. ([TOD and ( [Tl| ) respectively. The functions h^(kr) and h^(kr) 
are the Riccati-Hankel functions.c3 The boundary conditions that xi{'^'^E) is to satisfy are 
given in Eq. (|T2|) . If we write ai = 2iaii, these boundary conditions can be written in terms 
of the coefficients as 

ai2hj{Qa) + Pi2hY{Qa) = aiji{ka) (A2a) 

ai2hl'iQa) + A2 = (A2b) 

Tiihi{kh) + Ti2hi{kh) = ai2hi{Qh) + Pi2K{Qb) (A2c) 

J'ah+\kh) + J'nk'ikh) = 9.^ai2hi\Qh) + Pi2K'{Qb)) , (A2d) 

where ji = Yii^t - K)^ hl\kb) = ^^^l^^^b, and so on. Then 

ji{ka)hY'{Qa) - ^jl{ka)hY{Qa) 

ai2(k) = aAk)- ^ (Ada) 

' ^hi{Qa)hT\Qa)-hi\Qa)hi{Qa) 

3i(ka)h^'(Qa) - ^fAkaVhi (Qa) 
'hT{Qa)hi\Qa)-hl\Qa)hi{Qa) 
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^ „ , [htmk'm - fhrmkm]diika)hT\Qa) - fjl{ka)hT{Qa)] 

J^ii(k) = — ai(k) 

[h+{kb)h;\kb) - hl\kb)h-{kb)][h+{Qa)h^'{Qa) - hl'{Qa)h;{Qa)] 

\himk\kb) - Qk'{Qb)k{kb)][Uka)ht'{Qa) - fjl{ka)ht{Qa)] 

H — X X X X aiik) 

[hi{kb)hl'{kb) - hi'{kb)h-{kb)][h-{Qa)hi'{Qa) - h^' {Qa)hi {Qa)] 

(A3c) 

^ „ , [himhVm - ^hrmhimmka)hi\Qa) - ^f,{ka)k{Qa)] 

J^Mk) = — X X X X X 7-^ X aiik) 

[h-{kb)h+\kb)-hl\kb)hi{kb)][hi{Qa)h-'{Qa)-hi'{Qa)hl{Qa)] 

ikmhrm - fhrmhimmka)hrm - ^f,{ka)k{Qa)] 

H — 7^ 7. «/ «; • 

[h-{kb)hi'{kb) - h;'{kb)hl{kb)][h^{Qa)hl'{Qa) - h^'{Qa)hl{Qa)] 

(A3d) 

The Ith radial ket reads up to a normalization factor as 

ai{k)ji{kr) < r < a 

Xi{r- E) = { ai2{k)hl{Qr) + Pi2{k)hY{Qr) a<r<b (A4) 
Tii{k)hj{kr) + Ti2{k)hl{kr) b < r < oo . 



2. S-Matrix Approach 

Now we calculate the /th partial S matrix and its poles. Because 

hf ~ e^Kkr-i^m ^ r -> oo , (A5) 

the function J^nk Ki^y be interpreted as an outgoing wave with amplitude J^n and Ti2k 
may be interpreted as an incoming wave with amplitude J^i2- Thus the expression for the 
/th partial S matrix in the energy representation is 

J^i2[k) 



The resonances are associated with the poles of the analytic continuation of Eq. ( |A6|) . The 
functions J^n and Juo are analytic, and therefore the poles of the Si matrix ( |A6| ) are the zeros 
of its denominator .E^ Then, the condition !Fi2{k) = provides the Ith. resonance energies, 

[kmk'm - ^k'mkmmka)k\Qa) - ^j'i{ka)k{Qa)] 
- [kmk'm - ^k'mkmmka)k\Qa) - ^j'i{ka)k{Qa)] = o. 

(A7) 
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3. Gamow Vector Approach 



As we did for / = 0, we shall prove that the S-matrix poles are the same as the complex 
eigenvalues obtained from the purely outgoing boundary condition. 
The Gamow vector in the position representation satisfies 



r H ;r — ^ + V{r 



2m r dr'^ 2mr'^ 
If we write (r, 9, (f)\zR, I, m) := lxi{r; ZR)Yi^rn 



2m dr"^ 
The general solution of 

Xiir; zr) 



+ 



2mr^ 



Zr, I, m) = ZR{r, 6, (j)\zR, I, m) 
then the radial part fulfills 
+ y{.r) ] xi{r; zr) = ZRXi{r; zr) . 



(A8) 



(A9) 



IS 



aiih^{kr) + Pnhi (kr) < r < a 
anhjiQr) + (3i2hl{Qr) a<r <b 
J^iih^{kr) + J^i2h'[{kr) b < r < oo , 



(AlO) 



where k is the complex wave number (^) and Q is given by Eq. (|32D. The boundary 
conditions that the Gamow vectors satisfy are equivalent to the case / = 0. Now h^' plays 
the role of outgoing wave and hj' that of the incoming wave, 

(Alia) 
(Allb) 
(Allc) 
(Alld) 
(Alle) 

' ^ oo . (Allf) 
If we define = 2iaii, Eq. ( JAIID can be written in terms of the coefficients of ^r) as 



Xi{0 


Zr) 


= 




Zr) 


= Xiia+;zR) 




Zr) 


= x'iia+;zR) 


Xi{b- 


Zr) 


= Xiib+;zR) 


X'lib- 


Zr) 


= x'i{b+;zR) 


Xiir 


Zr) 


~ h^{kr) , 



ai2hl{Qa) + Pnhi {Qa) 
ai2h^'{Qa) + j3i2hj'{Qa) 



aiiiika) 
^aij'iika) 



Tiihiikh) = ai2hUQb)+Pi2KiQb) 
Taht\kh) = ^{ai2h^'{Qb)+Pi2h;\Qb)) . 



(Al2a) 
(Al2b) 
(Al2c) 
(Al2d) 



In the matrix representation Eq. ( |A12| ) looks like 



/ -ji{ka) 
^3'i{ka) 


V 



hi{Qa) h^iQa) \ 
h+'{Qa) hY'{Qa) 
-h+ikb) h+{Qb) h^iQb) 

-hvm fhi'iQb) fh-'m I 



a/2 
VA2/ 



/0\ 






(A13) 



20 



This system has non-trivial solutions iff the determinant of the coefficients is equal to zero, 



-ji{ka) 





hiiQa) 


hi (Qa) 







h+'iQa) 


K'iQa) 





-h+{kb) 


htm 


Km 







fhfm 


ik'm 



0. 



(A14) 



Straightforward calculations then lead to 

[htmht'ikb) - ^ht'mht{kb)][Uka)K'iQa) - ^jl{ka)hY (Qa)] 

- IKmhfikb) - ^h-'mht{kb)]\ji{ka)ht\Qa) - ^j[{ka)ht{Qa)] = 0. 

(A15) 

which was found in the previous section using the S matrix 



Equation (|A15|) is just Eq. 
approach. 



4. Green Function Approach 



The /th radial Green function satisfies 

hH{i + i) 



2m dr"^ ~'~ 2mr^ 
It expression is given by 



+ V{r) - E Gi{r, r'; E) = -6{r - r') . 



r (r r'- - Xi{r<; E)^i{ry, E) 

Lri[r,r — — — — 

h? W[xu W 



(A16) 



(A17) 



The function xi is the solution of the time independent Schrodinger equation that van- 
ishes at the origin. Therefore, xi{f]E) is given by Eq. ([A4D . The function ijji also satisfies 
Eq. (^, but with the boundary conditions 



iJi{a- 



htikr) 



r — > oo . 



Then 



ipi{r; k) 



aii{k)ht{kr) + bii{k)hf{kr) < r < a 
ai2ik)htm + bi2ik)htiQr) a<r <b 
ht{kr) b < r < oo , 



(A18a) 
(Al8b) 
(Al8c) 
(AlSd) 
(Al8e) 



(A19) 



where the a and b coefficients express the continuity conditions in Eq. ( |A18|) . The /th Green 
function is given by 
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The poles of Eq. ( |A20D are the zeros of its denominator .0 This is the same condition as was 
found in the S'-matrix approach and leads to the same result as the Gamow vector approach. 
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FIG. 1. The resonance energies of the square barrier potential. 
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FIG. 2. The resonance wave numbers of the square barrier potential. 
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FIG. 3. Deformation of the path of integration into the second sheet of the energy Riemann 
surface; (a) for the decaying states and (b) for the growing states. 
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